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Abstract. The static electric dipole polarizability of NaN clusters with even N has been calculated in a 
collective, axially averaged and a three-dimensional, finite-field approach for 2 < N < 20, including the 
ionic structure of the clusters. The validity of a collective model for the static response of small systems 
is demonstrated. Our density functional calculations verify the trends and fine structure seen in a recent 
experiment. A pseudopotential that reproduces the experimental bulk bond length and atomic energy levels 
leads to a substantial increase in the calculated polarizabilities, in better agreement with experiment. We 
relate remaining differences in the magnitude of the theoretical and experimental polarizabilities to the 
finite temperature present in the experiments. 

PACS. 36.40.-c Atomic and molecular clusters - 31.15.E Density functional theory in atomic and molecular 
physics - 33.15.Kr Properties of molecules, electric polarizability 



1 Introduction 

The measurement of the static electric polarizability of 
sodium clusters jl]] and its interpretation in terms of the 
jellium model [§] was one of the triggers for the research 
activities that today form the field of modern metal clus- 
ter physics. The first theoretical studies were followed by 
several others with different methods and aims: density 
functional calculations using pseudopotentials J|,|| or tak- 
ing all electrons into account |^] aimed at a quantitative 
description of the experimentally observed effects, semi- 
classical approaches || focused on size-dependent trends, 
and the static electric polarizability served to test and 
compare theoretical concepts 0,H,|). Recently, the field 
received new inspiration from a second experimental de- 
termination of the static polarizability of small, uncharged 
Na clusters [jl0|. 

Whereas a qualitative understanding of the experi- 
ments can be obtained with relatively simple models, a 
quantitative theoretical determination of the polarizabil- 
ity requires knowledge of the ionic and electronic config- 
urations of the clusters. Great effort has been devoted in 
the past to determine these [||[ll],[l^,|l3|. However, taking 
all ionic and electronic degrees of freedom into account 
in a three-dimensional calculation is a task of consider- 
able complexity. Therefore, most of these studies were re- 
stricted to clusters with not more than nine atoms. To 
reduce the computational expense, approximations for in- 
cluding ionic effects were developed l^,^5|,^6|, pj|p^ ,^9t. 
A second problem, however, is the great number of close- 
lying isomers that are found in sodium clusters. This effect 
is especially pronounced when stabilization of an overall 



shape through electronic shell effects is weak, i.e. for the 
"soft" clusters that are found between filled shells. In the 
present work we present calculations for the static electric 
polarizability that include the ionic structure in a realistic 
way. We take into account a great number of isomers for 
clusters with up to 20 atoms, especially for the soft clusters 
that fill the second electronic shell. The theoretical con- 
cepts that we used in this study are introduced in section 
||, where we also discuss the relevant cluster structures. In 
section ^| we present our results and compare with other 
calculations and experimental work. Our conclusions are 
summarized in section ^. 



2 Theoretical concepts 

The starting point for the theoretical determination of the 
polarizability of a cluster is the calculation of the ionic and 
electronic configuration of the ground-state and close ly- 
ing isomers. In the present work, this was done in two 
steps. First, we calculated low-energy ionic geometries for 
a wide range of cluster sizes with an improved version [ fl9| 
of the "Cylindrically Averaged Pseudopotential Scheme" 
(CAPS) ||. In CAPS the ions are treated fully three- 
dimcnsionally, but the valence electrons are restricted to 
axial symmetry. The cluster ground state is found by si- 
multaneously minimizing the energy functional with re- 
spect to the set of ionic positions (simulated annealing) 
and the valence-electron density. For the exchange and 
correlation energy we used the local-density approxima- 
tion (LDA) functional of Perdew and Wang (20|, and for 
the pseudopotential we employed the recently developed 
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phenomenological smooth-core potential that reproduces 
low temperature bulk and atomic properties fl9|| . Detailed 
comparisons with ab initio calculations have shown |l9f 
that CAPS predicts ionic geometries of sodium clusters 
rather accurately since truly triaxial deformations are rare. 
Furthermore, in a second step we performed fully three- 
dimensional (3D) Kohn-Sham (KS) calculations to check 
the ordering of isomers and to calculate polarizabilities 
without axial restriction on the electrons, and also in- 
cluded configurations from 3D geometry optimizations into 
our analysis as discussed below. 

Fig. ^ schematically depicts the most important ionic 
geometries for neutral clusters with even electron num- 
bers between 2 and 20. (We have calculated the polariz- 
abilities also for many further and higher isomers which, 
however, are not shown in Fig. [I] for the sake of clarity. 
They were omitted from the discussion since they do not 
lead to qualitatively different results.) For the small clus- 




Na16a Na16b 






Na18c 




Na 20 a Na20 b Na20 c 

Fig. 1. Cluster structures Na4 to Na2o- See text for discussion. 



ters Na2, Na4 and Nelq. many other theoretical predictions 
have been made |^,^.^0[|l^,^], and our geometries are in 
perfect agreement with them. In addition, due to the con- 
struction of the pseudopotential Jl9| , the bond lengths are 
close to the experimental ones, as e.g. seen in the dimer, 
where our calculated bond length is 5.78ao and the exper- 
imental one |2l| is 5.82ao- For Nas and Naio, our results 
are in agreement with 3D density functional calculations 
H§,[n|. For Nai2, we do not now of any ab initio calcu- 
lations. Therefore, besides two low-energy configurations 
from CAPS [(a) and (b)], we also included a locally re- 
optimized low-energy geometry from a 3D, Hiickel model 
calculation jL{| in our analysis (c). Our 3D calculations 
confirm our CAPS results and find structures (a) and (b) 
quasi degenerate with a difference in total energy of 0.05 
eV, whereas structure (c) is higher by 0.4 eV. Two of the 
three geometries considered for Nai4 [(b) and (c)] were 
also found very similar in 3D Hiickel model calculations 
(lIHI, and both CAPS and the 3D KS calculations find 
all of them very close in energy. For Nai6, we find as the 
CAPS-ground state structure (a), and in our 3D calcula- 
tions structure (b) is quasi degenerate with (a), whereas 
structures (c) and (d) are higher by 0.08 eV and 0.5 eV. 
Due to their very different overall shapes, these isomers 
span a range of what can be expected for the polarizabil- 
ity. For Nais and Na2o, our structures are again in close 
agreement with the 3D density functional calculation of 
p3[ , and all three structures are quasi degenerate. 

The static electric polarizability was calculated in two 
different ways. The first is based on a collective description 
of electronic excitations. It uses the well known equality 



a = 2m_i, 
which relates the negative first moment 



(1) 



m_i 
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of the strength function 



Sq(E) = l<HQ|0>| 2 <5(£„ - E - E), 



(2) 



(3) 



to the static electric polarizability a in the direction spec- 
ified by the external (dipole) excitation operator Q. 

In the evaluation of the strength function, the excited 
states \v) are identified with collective excitations. A dis- 
cussion of this approach can be found in [g2| . The collec- 
tive calculations were carried out using the cylindrically 
averaged densities and the "clamped nuclei approxima- 
tion" Q] , i.e. the ionic positions were taken to be the same 
with and without the dipole field. We have also checked 
this widely used approximation in the context of our stud- 
ies and find it well justified, as discussed below. 

The static polarizability can also be calculated directly 
from the derivative of the induced dipole moment fi in the 
presence of an external electric dipole field F ( "finite field 
method"): 



2Fi 



,J = x,y,z, 



(4) 
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where 



f i j (F) = -e / r 3 n(r,F)d 3 r + eZ^R 
•* R 



(5) 



for ions with valence Z . Here one has to make sure that the 
numerically applied finite dipole field F is small enough 
to be in the regime of linear response, but that it is on 
the other hand large enough to give a numerically stable 
signal. We have carefully checked this and found that the 
used field strengths between 0.00001e/ao 2 and 0.0005e/ao 2 
meet both requirements. Applied to the axial calculations, 
this approach allows to obtain the polarizability in the z- 
direction. By employing this method with the 3D KS cal- 
culations we have checked the influences of the axial av- 
eraging and the collective model on the polarizability and 
found that the z-polarizabilities from the axial and the 
3D finite- field calculations agree within 1% on the aver- 
age for the low-energy isomers. This shows that the axial 
averaging is a good approximation for the clusters dis- 
cussed here. The performance of t he collective model will 
be discussed below in Section 3.1. The orientation of our 
coordinate system was chosen such that the z-axis is in 
that principal direction of the tensor of inertia in which it 
deviates most from its average value. The average static 
electric polarizability 
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length (together with the compressibility and the atomic 
3s-level) when used with the LDA, and correspondingly 
results in a higher polarizability, in better agreement with 
experiment. It is further interesting to note that also the 
polarizabilities calculated with the empirical Bardsley po- 
tential ||, which was constructed to reproduce atomic 
energy levels, are noticeably higher than the BHS-based 
values. This shows that the cluster polarizability is also 
sensitive to atomic energy levels, and the fact that our 
values are closest to the experiment thus is a natural con- 
sequence of the combination of correct atomic energy lev- 
els and bond lengths. 



Table 1. Averaged static electric polarizability of small 
sodium clusters in A 3 . Brd: density functional (DF) calculation 
with empirical non-local Bardsley pseudopotential BHS: 
DF calculation with ab initio non-local Bachelet, Hamann, 
Schliiter pseudopotential Q. All el.: all-electron DF calcula- 
tion including gradient corrections to the exchange-correlation 
functional ^j. TM: DF calculation with Troullier-Martins non- 
local pseudopotential GAUSS.: DF calculation based on 
the GAUSSIAN94 program with SU basis set 0. Present: 
present work, values from three-dimensional approach. ExpR: 
recent experiment [JL0|. 



of course is independent of the choice of coordinate system. 
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3 Results 

3.1 Comparison of different theoretical results 

Since all density functional calculations that we know of 
agree on the geometry of the smallest sodium clusters, 
these clusters can serve as test cases to compare differ- 
ent theoretical approaches. In Table [j] we have listed the 
averaged static dipole polarizability as obtained in dif- 
ferent calculations, together with the value obtained in 
the recent experiment of Rayane et al. |l(J. All calcula- 
tions reproduce the experimental trend and give the cor- 
rect overall magnitude. But also, all calculations under- 
estimate the polarizability. The magnitude of this under- 
estimation, however, varies considerably for the different 
approaches. Whereas our results are closest to the exper- 
iment and close to the theoretical ones of Ref. ifTof , with 
the largest difference to the experiment being 8 % for Nas, 
a difference of 27 % is found for this cluster in the calcu- 
lation based on the ab initio Bachelet, Hamann, Schliiter 
(BHS) Pseudopotential A good part of this difference 
can be explained by comparing the bond lenghts of the 
clusters. The BHS pseudopotential considerably underes- 
timates the bond lengths leading to a higher electron 
density and a lower polarizability. Our empirical smooth- 
core pseudopotential, on the other hand, was constructed 
to reproduce the experimental low-temperature bulk bond 



From comparison with the calculations that went be- 
yond the LDA H,|l(J, it however becomes clear that the 
empirical pseudopotentials by construction " compensate" 
some of the errors that are a consequence of the use of the 
LDA. Therefore, it would be dangerous to argue that the 
inclusion of gradient corrections, which have been shown 
to increase the polarizability, could bring our calculated 
values in agreement with experiment: going beyond the 
LDA but keeping the empirical LDA pseudopotentials could 
lead to a double counting of effects. We therefore conclude 
that, on the one hand, a considerable part of the earlier 
observed differences between theoretical and experimen- 
tal polarizabilities can be attributed to effects associated 
with errors in the bond lengths or atomic energy levels, 
but on the other hand, further effects must contribute to 
the underestimation with respect to experiment. We will 
come back to this second point below. 

In Table ^ we have listed the polarizabilities of Na2 to 
Na2o for the geometries shown in Fig. [j} The left half gives 
the polarizability as computed from the 3D electron den- 
sity with the finite-field method, and the right half lists 
the values obtained in the axially averaged collective ap- 
proach. For the clusters up to Na 8 , the two methods agree 
well and the differences for the averaged polarizabilities 
are less than 1% for Na2 and Na§, and 3% for Na4 and 
Nag. This shows that the collective description is rather 
accurate, which is remarkable if one recalls that we are 
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dealing with only very few electrons. Beyond Nas, the dif- 
ferences are 6 % on the average, which is still fair, but ob- 
viously higher. This looks counter-intuitive at first sight, 
because the collective description should become better 
for larger systems. However, for N > 8 there comes an 
increasing number ofparticlc-hole states close to the Mie 
plasmon resonance p2| , leading to increasing fragmenta- 
tion of the collective strength. m_i and thus a is sensitive 
to energetically low-lying excitations since their energies 
enter in the denominator in Eq. ^|, and this can lead to an 
underestimation of the polarizability. 



A3 
for the cluster 

geometries of Fig. [I]. Left half: three-dimensional, finite field 
calculation. Right half: cylindrically averaged, collective calcu- 
lation. 



3.2 Comparison with experiments 

Fig. U shows a for our ground state structures as obtained 
in the axial, collective approach and the 3D finite-field 
calculations, in comparison to the two available sets of ex- 
perimental data. The absolute values for the experiments 
were calculated from the measured relative values with 
an atomic polarizability of 23. 6A jn|. To guide the eye, 
the polarizabilities from each set of data are connected 
by lines. Both experiments and the theoretical data show 
that, overall, the polarizability increases with increasing 
cluster size. The polarizability from the axial collective 
model qualitatively shows the same behavior as the one 
from the 3D finite field calculation. Comparison of the 
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29.9 
47.0 
129.8 
118.1 
126.3 
213.4 
156.5 
158.6 
183.5 
271.3 
175.2 
193.5 
235.2 
212.5 
272.2 
260.3 
251.8 
281.6 
285.7 
275.0 
267.9 



29.9 
59.0 
129.8 
118.5 
126.3 
213.4 
158.6 
145.1 
183.5 
273.9 
179.3 
193.5 
235.2 
213.4 
260.8 
262.5 
250.7 
280.5 
284.5 
275.0 
271.5 



51.3 
130.1 

62.2 
132.3 
219.3 
155.7 
261.1 
285.3 
278.9 
137.8 
291.3 
394.5 
231.5 
272.2 
239.3 
291.3 
283.8 
250.7 
309.4 
311.8 
295.2 



37.0 
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239.3 
271.4 
262.1 
270.9 
293.2 
287.3 
278.2 



30.2 
53.7 
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117.5 
125.3 
199.5 
151.3 
148.7 
177.5 
264.2 
171.1 
190.3 
227.3 
206.0 
260.6 
236.4 
232.6 
255.9 
269.7 
261.8 
283.2 



51.5 
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62.1 
131.4 
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132.7 
262.5 
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255.1 
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266.7 
253.4 
227.7 
279.7 
282.9 
280.3 
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76 
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148 
180 
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220 
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233 
225 
222 
235 
246 
239 
246 
273 
268 
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Fig. 2. Static electric dipole polarizability in A versus number 
t oi electrons. Crosses with thin, long dashed line: experiment of 
^Rayane et al. [ju^; stars with strong, long dashed line: exper- 
iment of Knight et al. Q; open squares with full line: present 
^work, three-dimensional finite-field calculation for lowest iso- 
'mer; filled squares with short dashed line: present work, axially 
''averaged collective calculation for lowest isomer. See text for 
^discussion of error bars. 



Comparing the polarizabilities of clusters with the same 
number of electrons but different geometries shows the in- 
fluence of the overall shape of the cluster. For Nai4, e.g., 
isomers (a) and (c) have a valence electron density which 
is close to prolate, whereas (b) has a more oblate one. The 
averaged polarizability for the two prolate isomers is equal, 
although their ionic geometries differ. The oblate isomer, 
however, has a noticeably higher averaged polarizability. 
This is what one expects, because for oblate clusters there 
are two principal directions with a low and one with a high 
polarizability, whereas for prolate clusters the reverse is 
true. The fact that different ionic geometries can lead to 
very similar averaged polarizabilities is also seen for Nai2- 
It thus becomes clear that contrary to what was believed 
earlier Q one cannot necessarily distinguish between de- 
tails in the ionic configuration by comparing theoretical 
values to experimental data that measure the averaged 
polarizability. 



3D values with the experimental data shows that for the 
smallest clusters, the theoretical and experimental values 
agree as discussed before, and the values obtained in the 
two experiments are comparable up to Naio ■ Beyond Naio i 
the discrepancies between the two experiments become 
larger, and also the differences between theoretical and ex- 
perimental polarizabilities increase. For Nai2, Nai4, Nai6 
and Nai8 the experiment of Knight et al. gives lower values 
than the experiment of Rayane et al., and the calculated 
averaged polarizability is lower than both experiments for 
Nai2, Nai4, and Nai 6 . For Nai 8 the finite-field value ob- 
tained for our ground-state structure matches the value 
measured by Knight el al., and for Na2o, our ground-state 
polarizability is very close to the measurement of Rayane 
et al. In this discussion one must keep in mind, however, 

that the experimental uncertainty is about +/-2 A per 
atom ]l(]|] , i.e. the uncertainty in the absolute value in- 
creases with the cluster size, as indicated by the errors 
bars in Fig. El Comparisons are made easier if the linear 
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Fig. 3. Normalized static electric dipole polarizability. Crosses 
with thin dashed line: experiment of Rayane et al. jtol; stars 
with strong dashed line: experiment of Knight et al. llJTsquares 
with full line: present work, three-dimensional finite-field cal- 
culation for lowest isomer; open triangle, filled triangle and 
upside-down triangle: second, third and fourth isomer, respec- 
tively; filled circles: Jellium results from J2,R| . See text for dis- 



growth in a is scaled away. Therefore, one should rather 
look at the normalized polarizability 



Na a 



(7) 



which is shown in Fig. ||, because it allows to identify 
trends and details more clearly. 

From Fig. || it becomes clear that for Na2 to Nas, the 
trend seen in the two experiments is similar up to one 
exception: For Na@, the experiment of Rayane et al. pre- 
dicts a noticeably smaller value than the one by Knight 
et al. Comparison with our theoretical data shows that, 
although the values of the older experiment are closer to 
the theory with respect to magnitude for Na2, Na4 and 
Nas, the trend that our data show corresponds clearly to 
the one seen in the new experiment since the two curves 
are parallel. Going from Nas to Naio, both experiments 
predict a steep rise in the polarizability. This rise due to 
the shell closing at Nas is also seen in the theoretical data, 
but it is less pronounced than in the experiments (as we 
will discuss below). For Nai2, a higher a 11 than for Naio 
is predicted by the data of Rayane, whereas the reverse 
ordering is seen in the data of Knight et al. Again, our 
calculations support the finding of the new experiment, 
and all isomers lead to similar a n . For Nai4, both experi- 
ments show a decrease. Our prolate ground state and iso- 
mer reproduce this trend. That it is the prolate structures 
that fit to the experiment is consistent with the ah initio 
molecular dynamics calculations of Hakkinen et al. f23f . 
The next step to Nai6 again reveals a slight difference 
between the two experiments: both predict an increase 
compared to Nai4, but whereas the older experiment sees 
a n smaller for Nai6 than for Nai2, the new experiment 
shows the opposite ordering. Once more, our ground state 



structure leads to a polarizability that follows the trend of 
the new experiment. (The other isomers, however, lead to 
smaller polarizabilities, and an explanation for the differ- 
ence between the experiments thus might be that different 
ensembles of isomers were populated due to slightly dif- 
ferent experimental conditions.) Going to Nais leads to 
a decrease in the polarizability in both experiments. Our 
calculation shows this decrease, which is a manifestation of 
the nearby shell closing. But whereas the old experiment 
actually sees the shell closing at Nais and an increase in 
the polarizability for Na2o, the new experiment and our 
data find an absolute minimum at Na2o- 

A comparison with the polarizability obtained in the 
spherical jellium model ||,|| for Nas and Na2o, also indi- 
cated in Fig. ^, shows the improvement that is brought 
about by the inclusion of the ionic structure. 

3.3 Discussion 

As just discussed, our calculations reproduce the fine struc- 
ture seen in the new experiment. However, there is no 
obvious explanation for why the results of the two exper- 
iments differ p4| . Also, there is a characteristic change 
in the magnitude of the difference between our theoreti- 
cal results and the experimental values of Rayane et al.: 
whereas the calculated values for Na2 to Nas and Na2o on 
the average differ only by 5 % from the new experiment, 
the open-shell clusters from Naio to Nais show 18 % dif- 
ference for the ground state. The increase in polarizability 
when going from Nas to Naio is considerably underesti- 
mated, whereas the following steps in the normalized po- 
larizability are nearly reproduced correctly, i.e. it looks as 
if the theoretical curve for Naio to Nais should be shifted 
upwards by a constant. A first suspicion might be that this 
"offset" could be due to the use of CAPS in the geometry 
optimization. But it should be noted that the step occurs 
at Naio, and that also Nai4 and Nais are off by the same 
amount. Since the low-energy structures of these clusters 
are well established, as discussed in Section |[ and since 
also the geometry for Nai2 from the 3D calculation does 
not lead to qualitatively different results, we can conclude 
that the differences are not due to limits in the geometry 
optimization with CAPS. Also, the neglection of the re- 
laxation of the nuclei in the presence of the electric dipole 
field has been investigated earlier j| for the small clus- 
ters and was shown to be a good approximation. We have 
counter checked this result for the test case Naio and find 
corrections of less than 1%. 

An obvious limitation of our approach is the neglect 
of the core polarization. However, the all-electron calcu- 
lations of Guan et al. || treat the core electrons explic- 
itly and do not lead to b ette r agreement with experiment, 
as discussed in Section 3.1. From this, one already can 
conclude that core polarization cannot account for all of 
the observed differences. Its effect can be estimated from 
the polarizability of the sodium cation. Different measure- 
ments find values between 0.179 A and 0.41 A , lead- 
ing to corrections of - roughly - 1-2 % in a n . Since the core 
polarizability leads to a shift in a n that is the same for all 
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cluster sizes, it contributes to the difference that is also 
seen for the smallest clusters, but it cannot explain the 
jump in the difference seen at Naio- 

Another principal limitation of our approach is the use 
of the LDA. As, e.g., discussed in ||, the LDA can af- 
fect the polarizability in different and opposing ways. On 
the one hand it may lead to an overscreening and thus 
an underestimation of the polarizability, and early calcu- 
lations within the spherical jellium model reported that 
indeed the static polarizability was increased if one went 
beyond LDA using self- interaction corrections (tJ . On the 
other hand, self-interaction corrections can lead to more 
negative single particle energies and thus to smaller polar- 
izabilities [||, and Refs. ||] and |2(| give examples where 
the overall effect of self-interaction corrections on the op- 
tic response is very small. One cannot directly conclude 
from the jellium results to our ionic structure calculations, 
because the sharp edge of the steep- wall jellium model can 
qualitatively lead to differences. But in any case it is highly 
implausible that the LDA affects the clusters from Naio 
to Nais much stronger than the other ones, and it should 
also be kept in mind that the worst indirect effects of the 
LDA as, e.g., underestimation of bond lengths, are com- 
pensated by using our phenomenological pseudopotential. 

One might also ponder about possible uncertainties in 
the experimental determination of the polarizabilities. A 
considerable underestimation could be explained if one as- 
sumes that while passing through the deflecting field, the 
clusters are oriented such that one always measures the 
highest component of the polarizability. In that case, we 
would not have to compare the averaged value to the ex- 
periment, but the highest one. One could imagine that 
the cluster's rotation be damped, since angular momen- 
tum conservation is broken by the external field and the 
energy thus could be transferred from the rotation to in- 
ternal degrees of freedom (vibrations). The time scale of 
this energy transfer is not known, but since the clusters 
are spending about 10 _4 s in the deflecting field region, it 
seems unlikely that there should be no coupling over such 
a long time. One could further argue that for statistical 
reasons it is less likely that a larger cluster will lose its 
orientation again through random-like thermal motion of 
its constituent ions than a smaller cluster. However, the 
maximal energy difference between different orientations 
is very small, for Naio, e.g., it is 0.3Kk_R for the typical 
field strength applied in the experiment Thus, thermal 
fluctuations can be expected to wipe out any orientation. 

From another point of view, however, the finite tem- 
perature explains a good part of the differences that our 
calculations (and other calculations for the small clusters) 
show in comparison to the experimental data. Whereas 
our calculations were done for T=0, the supersonic nozzle 
expansion used in the experiment produces clusters with 
an internal energy distribution corresponding to about 400 
- 600 K [^,|2^]. An estimate based on the thermal ex- 
pansion coefficient of bulk sodium leads to an increase 
in the bond lengths of about 3 %, and a detailed finite- 
temperature CAPS calculation |^9| for Na^ at 400 K also 
shows a bond length increase of 3 %. This will only be a 



lower limit, since in neutral clusters one can expect a larger 
expansion than in the bulk due to the large surface, and 
also a larger expansion than for charged clusters. Thus, to 
get an estimate for the lower limit of what can be expected 
from thermal expansion, we have scaled the cluster coor- 
dinates by 3 % and again calculated the polarizabilities, 
finding an increase of about 3 % for the planar and 5 % 
for three-dimensional structures. This finding is consistent 
with the results of Guan et al. Together with the correc- 
tions that are to be expected from the core polarizability, 
this brings our results for the small and the closed shell 
clusters in quantitative agreement with the experimental 
data. 



4 Summary and Conclusion 

We have presented calculations for the static electric dipole 
polarizability for sodium clusters with atom numbers be- 
tween 2 and 20, covering several low-energy structures 
for each cluster size beyond Naio. By comparing our re- 
sults to previous calculations for the smallest clusters, we 
have shown that a pseudopotential which correctly repro- 
duces atomic and bulk properties also improves the static 
response considerably. We have shown that a collective 
model for the excited states of sodium clusters, whose va- 
lidity for the dynamical response was established previ- 
ously, works reasonably also for the static response in real- 
istic systems. Over the whole range of cluster sizes studied 
in the present work, we confirm the fine structure seen in a 
recent experiment. By comparing the calculated averaged 
polarizability of different isomers for the same cluster size 
to the measured polarizability, we showed that completely 
different ionic geometries can lead to very similar averaged 
polarizabilities. By considering higher isomers we further- 
more took a first step to take into account the finite tem- 
perature present in the experiment. Our results show that 
for the open shell clusters from Naio to Naig, also higher 
lying isomers do not close the remaining gap between the- 
ory and experiment. This shows that it is a worthwhile 
task for future studies to investigate the influence of finite 
temperatures on these "soft" clusters explicitly. For Na2 
to Nag and Na2o, we showed that quantitative agreement 
is already obtained when the effects of thermal expansion 
and the core polarizability are taken into account. 



One of us (S. Kiimmel) thanks K. Hansen for several clar- 
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